Implementation of the quantum walk step operator in lateral quantum dots 
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We propose a physical implementation of the step operator of the discrete quantum walk for an 
electron in a one-dimensional chain of quantum dots. The operating principle of the step operator is 
based on locally enhanced Zeeman splitting and the role of the quantum coin is played by the spin 
of the electron. We calculate the probability of successful transfer of the electron in the presence of 
decoherence due to quantum charge fluctuations, modeled as a bosonic bath. We then analyze two 
mechanisms for creating locally enhanced Zeeman splitting based on, respectively, locally applied 
electric and magnetic fields and slanting magnetic fields. Our results imply that a success probability 
of > 90% is feasible under realistic experimental conditions. 
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The quantum walk (or quantum random walk) is the 
quantum-mechanical analogue of the classical random 
walk and describes the random walk behavior of a quan- 
tum particle. The concept "quantum walk" was formally 
introduced by Aharonov et al. in 1993 [l[ and suggested 
earlier by Feynman [2J. The essential difference with the 
classical random walk lies in the role of the coin: Whereas 
in the classical random walk the coin is a classical ob- 
ject with two possible measurement outcomes ("heads" 
or "tails"), in the quantum walk the coin is a quantum- 
mechanical object - typically a two-level system such as a 
spin- 1/2 particle - which can be measured along different 
bases and hence has a multi-sided character. As a result, 
the quantum walk exhibits strikingly different dynamic 
behavior compared to its classical counterpart due to in- 
terference between different possible paths. One exam- 
ple is faster propagation [l[: The root-mean-square dis- 
tance from the origin (x) Tms that is covered by a quantum 
walker grows linearly with the number of steps N (thus 
corresponding to ballistic propagation), whereas for the 
classical random walk (x) rms ~ \/~N (corresponding to 
diffusive propagation) . This property has been exploited 
to design new quantum computing algorithms Both 
discrete and continuous time quantum walks have been 
extensively studied in recent years including inves- 
tigations of decoherence [B| and entanglement between 
quantum walkers Q. 

As far as implementations of quantum walks in ac- 
tual physical systems are concerned, several proposals 
have been put forward for a range of optical and atomic 
systems, such as optical cavities ml, cavity QED sys- 
tems 0, trapped atoms and ions 9] and linear optical 
elements [l(J. On the experimental side, only a few re- 
alizations of quantum walks have been achieved: Dis- 
crete and continous quantum walks in NMR quantum 
systems II ill, discrete quantum walks using linear optical 
elements [12j | and, most recently, a continuous quantum 
walk in an optical waveguide lattice [3] • 



For solid-state systems, no realizations of quantum 
walks exist so far. A recent proposal for implementation 
of a quantum algorithm using NAND operations in a tree 
of quantum dots relies on the continuous time quantum 
walk 



In this paper we propose the first implementation of 
a discrete quantum walk in a solid-state quantum sys- 
tem, which consists of a single electron traveling in a one- 
dimensional chain of quantum dots. In particular, we fo- 
cus on the implementation of the so-called step operator, 
the basic unit of the quantum walk. The step operator 
causes the electron to either move to the left or to the 
right depending on the state of the quantum coin, which 
in our model is represented by the spin of the electron. 
We calculate the spin-dependent transfer probability of 
the electron from one dot to the next in the presence of 
different energy level splittings in neighboring quantum 
dots (due to locally enhanced Zeeman splitting), tak- 
ing into account the effects of decoherence due to gate 
voltage fluctuations [lj|. We then propose two physical 
mechanisms to achieve locally enhanced level splitting in 
a quantum dot using local electric and magnetic fields 
and find that under current experimental circumstances 
successful implementation of the step operator is possible 
with > 90% probability. 

Model of the step operator. - Consider a chain of three 
quantum dots in series, in which the middle dot (M) is 
occupied by a single electron, see Fig. Q] The energy 
level diagram shows the Zeeman splitting A{' ^ of the low- 
est two levels in a magnetic field, and we assume that 
this splitting is larger in the middle dot by an amount 
A' z — A z . The gate voltages Vl, Vm and Vr are used 
to shift the energy levels in the left, middle and right 
dot resp., and the gate voltages Vr are used to tune the 
tunnel coupling between neighboring dots. The initial 
spin state of the electron is a superposition of spin-t and 
spin-| and our first goal is to design an implementation 
of the step operator such that it causes the electron to 
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Magnetic and Eleti.ru' Fields 

FIG. 1: (left) A linear chain of 3 quantum dots in a magnetic 
field and (right) the corresponding energy level diagram. The 
direction of (locally) applied electric and magnetic fields in 
each dot is also indicated, a cross indicates no applied field. 
See the text for further explanation. 



coherently move to the left (right) if it has spin-]" (spin- J,) 
and to calculate the probability of succesful transfer for 
this process. To begin with, we assume the right dot to 
be decoupled from the other two and calculate the spin- 
dependent probabilities to find the electron in the left and 
middle dot at a given time. In the absence of decoherence 
- which is considered further below - the Hamiltonian for 
the spin-| and spin- J, components of the electron is given 
by 



H 



0,(T 



X 



X 



(1) 



in the basis {\L),\M)}. Here (A u B r ) = (e,0), 
= (e + A 2 ,A Z ) and A is the (spin-independent) tun- 
nel coupling between the left and middle dot. The 
eigenvalues and -vectors of Eq. (P) are given by 
E±, a = ± y(A a -B a y + AX\ |^ +)CT ) = 
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(sin# CT , cos 9 a) and \tp- 



(cos 9 a , — sin 9 a ) , with 

tan 9 a = v 2A , Of = e/Ti and 6^ = (e + A z — 

A z )/h. In the absence of decoherence, the solution of 
the density matrix equations p = —(i/h)[Ho^, p] for the 
population in the left dot is given by 



PLL,a{t) 
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[1 — cos(w CT t)] 
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for initial conditions pll,o-{0) = 0, pmm,<t(0) = 1 and 
Tiu! a = sj (hd a ) 2 + AX 2 . We see that the probability 
for the spin-| component to be in the left dot and the 
spin-| component to remain in the middle dot is 1 for 

= 1 - cos hmr^j , n = 0, 1, 2, . . .. The smallest 

n that yields a solution to this equation for e = (zero 
detuning) is n = 1, for which A z — A' z = 2^3 X and 
t = 27t/oj|. The second half of the step operator then 
consists of repeating the same procedure as described 
above for tunnel coupling between the middle and the 
right dot, now assuming the left dot to be decoupled. 

Decoherence due to quantum fluctuations. - In reality, 
the time evolution of the occupation probabilities p a is 
affected by decoherence due to coupling of the quantum 



dots to the environment. In typical experimental situa- 
tions kuT < hw„, for T < 100 mK and A < 1.5 10~ 24 
Il7j ]. so that quantum noise, rather than classical 
noise, is the dominant source of decoherence. Specifi- 
cally, since the quantum walk step operator based on the 
Hamiltonian JTJ) involves tuning the tunnel couplings A - 
resulting in different charge occupations on neighboring 
dots - the probability distributions will be strongly af- 
fected by fluctuations of charge in the environment [3] • 
The quantum charge fluctuations we consider here are 
gate voltage fluctuations, which cause both fluctuations 
in the tunnel coupling A and in the energy levels E± of 
the Hamiltonian H) 19]. The Hamilt onian which de- 
scribes the quantum dot system plus the environment is 
given by 



H = H + V e A e + V X A X 
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with Ho the Hamiltonian ([T]) of the isolated sys- 
tem, V x = \L)(M\ + \M)(L\, V e = \L)(L\ - \M)(M\, 

A e(X) = J2k C kA^)( b Le(\) + & M(A)) and #&atft,e(A) = 

J2k fy^k,e(x)^\ e (A)kfc,e(AV We model the environment 
as a bosonic bath 12011 with creation and annihila- 



tion operators 



6 f 

fc,e(A) 



and b 



fc,e(A) 



and use the Born- 



Markov approximation [21[ to calculate the time evolu- 
tion of the spin-dependent occupation probabilities un- 
der the Hamiltonian ([3|), assuming weak coupling be- 
tween the system and the environment and short cor- 
relation times of the boson baths (Markovian assump- 
tion). The baths are characterized by the symmet- 
ric and anti-symmetric spectral functions S (oj) with 

S + {uj) = coth (^t) ^"M and we assume Ohmic 
baths with Lorentzian damping for which S~(u>) = 
ah 2 uj ( w /J~yi + i with lo c a cut-off frequency. The time 
evolution of the reduced density matrix p a t(t) is then 

given by dp °f t = ~iuJ a bPab{t) + X^ s rf cc) RabcdPcd{t) , with 
uj ab = (E n — E b )/h and R a bcd the Bloch-Rcdficld ten- 
sor [2C|, 22] . Solving this master equation for the two-level 
system {|i/'+, CT ), \ip-, a )} described above for an electron 
starting at t = in the middle dot yields for the evolu- 
tion of the population p of the groundstate \ip~) and 

the coherence terms p |_: 
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FIG. 2: Survival probability [Eq. ©J of the spin-|-component 
(main plot) and spin-|-component (inset) in the middle dot 
as a function of time t for various values of a = a\ — 
a e Parameters used are h8^ = 0, hS^ = 1/ieV, A = 



(2>/3) /ieV [H|], and T = 10 mK. 



Eqns. m and ((5|) are valid in the limit 71, 72 <C 
I TioJu 1 , 1 /f (with f the bath correlation time) [21| . The 
survival probability to remain in the middle dot at time 
t is then given by 



1 + cos(2^o.) tanh 



( hu 



\2k B T 
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and plotted for both spin-]" and spin- j in Fig. [2] Using 
typical experimental parameters (see figure caption and 
discussion below) we see that succesful transfer of the 
spin-| component to the left dot while spin-J, remains 
in the middle dot occurs with probability > 90% (for 
a < 3 • 1CT 4 ) at * » 3.5 ns. 

Locally enhanced Zeeman splitting. - We now consider 
the question how a different Zeeman splitting in neigh- 
boring dots - as assumed in the calculations above - can 
be realized in practice. In particular, we propose two 
mechanisms to achieve locally enhanced Zeeman split- 
ting. The first one is by application of a local trans- 
verse magnetic field and a local electric field and relies 
on spin-orbit interaction. Both tunable local magnetic 
and electric fields have recently been demonstrated ex- 
perimentally 2J, |25[ . In our model, see Fig. [TJ we as- 
sume that a local magnetic field (in addition to the global 
Zeeman-splitting field B) and a local electric field E(t) 
are applied to the middle quantum dot. In the presence 
of spin-orbit interaction, the Hamiltonian for this middle 
dot is then given by H = Ho+eE(t)-r+a (p x a y — p y (T x ) + 
(3 (— p x a x + pyUy), with Hq the Hamiltonian given by 
Eq. ([TJ) , and a and (3 resp. the Rashba and Dressel- 
haus spin-orbit coupling strengths. Using the Schrieffer- 
Wolff transformation, the Hamiltonian H can be di- 
agonalized to first order in the spin-orbit terms which 
yields, for the spin-dependent part, H c g = e s He~ s — 
^9*f l B (b + SB(t)j ■ a [26] . Here the effective magnetic 

field SB is given by SB = 2B x (fii(i) +11X n 2 (tf 



with f2i(i) 



§-oci{E y ,/\_,E x ,/\ + ,Q), n 2 (t) 
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Ih 
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3r, A± = h/(m*((3±a)), E z 



g*HsB, E ia = Ti{yjAujl + w 2 ±w c )/2, with lu c = eB/m* 
the cyclotron frequency and m*Lo^r 2 jl the harmonic po- 
tential of the quantum dot In these expressions we 
have used x' = (x + y)/ \[2 and y' = (y — x)/ \/2. Since 
for the step operator we require the eigenstates of the 
Hamiltonian H to be | f) and | J.) we must account for 
the fact that the eigenvectors of the effective Hamiltonian 
H e s are transformed by e s with respect to the eigenvec- 
tors of H corresponding to the same eigenvalues. We do 
this be requiring that the eigenvectors of H c g are e s \ f) 
and e s \ J.) to first order in the action S. 

Let us assume the global magnetic field B to be paral- 
lel to the z-axis. From the expression for H c s it follows 
that in order to locally generate a different Zeeman split- 
ting (along the z-axis) in the middle dot compared to its 
neighboring dots, an additional tunable local magnetic 
field, e.g. along the i'-axis, and a tunable local electric 
field, also along the i'-axis, are required. We define the 
total (local) magnetic field B(t) — B (n x '(t),0,n z ). We 
use the expression for the effective magnetic field com- 
bined with the requirement for the eigenvectors of H c h 
to derive two implicit equations for the magnetic- and 
electric field in the i'-direction as a function of the local 
magnetic field SB in the z-direction: 
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where we have defined 7 = 2eh~/Ez and 
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Fig. [3] shows the required fields B x > and E x > as a func- 
tion of Bo for typical experimental parameters and a local 
Zeeman splitting of 1 \ieV [28| . We see that in order to 
have an additional Zeeman splitting of 1 /j,eV (as assumed 
in Fig. [5]), e.g. at B = 0.05T magnetic and electric fields 
B x > ~ 63 mT and E x > ~ 4.6 10 5 Vm" 1 are required. 
Although these values of B x i and E x > are (somewhat) 
larger than those that have been used so far in exper- 
iments [3, l30j . they may well become available in the 
near future. Alternatively, one could use larger Bq (cor- 
responding to larger B x i and smaller E x i ) or smaller SB Z 
(corresponding to smaller B x i and E x i , but also smaller 
success probability). In addition, we note that in order 
to preserve the qubit's quantum state, the fields need to 
be switched on and off adiabatically. The correspond- 
ing switching time T of e.g. the tunable magnetic field 



then has to fulfill 311 T > 



0.1 ns. This 



in^\g'\ii B B 

timescale is well within experimental reach and compat- 
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dots. For currently available techniques to locally create 
enhanced Zeeman splitting and taking into account de- 
coherence due to quantum charge fluctuations, we have 
analyzed the probability for coherent transfer of the elec- 
tron to the left or right (conditioned on its spin), and pre- 
dict that > 90 % success probability is feasible. We hope 
that our results will stimulate a proof-of-principle exper- 
imental demonstration of the discrete quantum walk step 
operator in a solid-state nanosystem. 

This work has been supported by the Netherlands Or- 
ganisation for Scientific Research (NWO). 



FIG. 3: E x i (solid line) and B x i (dotted line) vs. Bo for a 
local magnetic field 8B corresponding to a Zeeman splitting of 
1 fieV. Parameters used are Tiloq = 1 me V |27| . m* = 0.067m e 
(for GaAs [I]}), A+ = 2^m, A_ = 10 /im [MlS], and n z = 1. 



ible with the operation time (duration of electron trans- 
fer) w 3.5 ns that we estimated above. 

Another method to generate different Zeeman fields in 
neighboring dots is by using a slanting magnetic field. 
The latter has recently been demonstrated for the first 
time by integrating a microsize ferromagnet in a dou- 
ble dot device [32]. As a result of the slanting field, the 
orbital and spin degrees of freedom become hybridized, 
leading to an effective mixed charge-spin two-level sys- 
tem, where the role of the coin is played by the pseu- 
dospin instead of the real spin. A global magnetic field 
~ 2 T magnetizes the ferromagnet and its inhomogene- 
ity leads to a different Zeeman field in the two quantum 
dots of SB Z ~ 10 mT (|A 2 - A' z \ w 4- 10~ 26 J [32]). An 
advantage of this method compared to using spin-orbit 
interaction (as discussed above) is that no tunable mag- 
netic fields but only tunable gate voltages are needed. A 
disadvantage is that the qubit is likely to be more sensi- 
tive to orbital decoherence. 

Finally, we briefly discuss the question of how to ex- 
tend our model from 3 quantum dots to a longer one- 
dimensional chain. In order for the electron to perform 
a quantum walk along the chain, opening and closing 
of tunnel barriers and aligning energy levels has to be 
applied at every position where the particle has a finite 
probability of being found. In addition, not only a step 
operator, but also a reinitialization operator C of the 
coin (spin) degree of freedom is needed EL A common 
choice for C is the Hadamard operator [J], which can 
be implemented by two coherent rotations around differ- 
ent axes [iil I25L l33l | . The dynamics of a quantum walk 
of electrons along a one-dimensional chain of quantum 
dots in the presence of decoherence remains an interest- 
ing question for future research. 

In conclusion, we have proposed an implementation of 
a discrete quantum walk step operator in a solid-state 
nanostructure consisting of a linear chain of 3 quantum 
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